Changes in environmental exposures over decades may influence the genetic architecture of severe spermatogenic failure.

STUDY QUESTION: Do the genetic determinants of idiopathic severe spermatogenic failure (SPGF) differ between generations? SUMMARY ANSWER: Our data support that the genetic component of idiopathic SPGF is impacted by dynamic changes in environmental exposures over decades. WHAT IS KNOWN ALREADY: The idiopathic form of SPGF has a multifactorial etiology wherein an interaction between genetic, epigenetic, and environmental factors leads to the disease onset and progression. At the genetic level, genome-wide association studies (GWASs) allow the analysis of millions of genetic variants across the genome in a hypothesis-free manner, as a valuable tool for identifying susceptibility risk loci. However, little is known about the specific role of non-genetic factors and their influence on the genetic determinants in this type of conditions. STUDY DESIGN, SIZE, DURATION: Case-control genetic association analyses were performed including a total of 912 SPGF cases and 1360 unaffected controls. PARTICIPANTS/MATERIALS, SETTING, METHODS: All participants had European ancestry (Iberian and German). SPGF cases were diagnosed during the last decade either with idiopathic non-obstructive azoospermia (n = 547) or with idiopathic non-obstructive oligozoospermia (n = 365). Case-control genetic association analyses were performed by logistic regression models considering the generation as a covariate and by in silico functional characterization of the susceptibility genomic regions. MAIN RESULTS AND THE ROLE OF CHANCE: This analysis revealed 13 novel genetic association signals with SPGF, with eight of them being independent. The observed associations were mostly explained by the interaction between each lead variant and the age-group. Additionally, we established links between these loci and diverse non-genetic factors, such as toxic or dietary habits, respiratory disorders, and autoimmune diseases, which might potentially influence the genetic architecture of idiopathic SPGF. LARGE SCALE DATA: GWAS data are available from the authors upon reasonable request. LIMITATIONS, REASONS FOR CAUTION: Additional independent studies involving large cohorts in ethnically diverse populations are warranted to confirm our findings. WIDER IMPLICATIONS OF THE FINDINGS: Overall, this study proposes an innovative strategy to achieve a more precise understanding of conditions such as SPGF by considering the interactions between a variable exposome through different generations and genetic predisposition to complex diseases. STUDY FUNDING/COMPETING INTEREST(S): This work was supported by the "Plan Andaluz de Investigación, Desarrollo e Innovación (PAIDI 2020)" (ref. PY20_00212, P20_00583), the Spanish Ministry of Economy and Competitiveness through the Spanish National Plan for Scientific and Technical Research and Innovation (ref. PID2020-120157RB-I00 funded by MCIN/ AEI/10.13039/501100011033), and the 'Proyectos I+D+i del Programa Operativo FEDER 2020' (ref. B-CTS-584-UGR20). ToxOmics-Centre for Toxicogenomics and Human Health, Genetics, Oncology and Human Toxicology, is also partially supported by the Portuguese Foundation for Science and Technology (Projects: UIDB/00009/2020; UIDP/00009/2020). The authors declare no competing interests. TRIAL REGISTRATION NUMBER: N/A.

Non-classical circulating monocytes expressing high levels of microsomal prostaglandin E2 synthase-1 tag an aberrant IFN-response in systemic sclerosis.

Systemic sclerosis (SSc) is a complex disease that affects the connective tissue, causing fibrosis. SSc patients show altered immune cell composition and activation in the peripheral blood (PB). PB monocytes (Mos) are recruited into tissues where they differentiate into macrophages, which are directly involved in fibrosis. To understand the role of CD14+ PB Mos in SSc, a single-cell transcriptome analysis (scRNA-seq) was conducted on 8 SSc patients and 8 controls. Using unsupervised clustering methods, CD14+ cells were assigned to 11 clusters, which added granularity to the known monocyte subsets: classical (cMos), intermediate (iMos) and non-classical Mos (ncMos) or type 2 dendritic cells. NcMos were significantly overrepresented in SSc patients and showed an active IFN-signature and increased expression levels of PTGES, in addition to monocyte motility and adhesion markers. We identified a SSc-related cluster of IRF7+ STAT1+ iMos with an aberrant IFN-response. Finally, a depletion of M2 polarised cMos in SSc was observed. Our results highlighted the potential of PB Mos as biomarkers for SSc and provided new possibilities for putative drug targets for modulating the innate immune response in SSc.

Celiac Disease Is a Risk Factor for Mature T and NK Cell Lymphoma: A Mendelian Randomization Study.

Celiac disease (CeD) is an immune-mediated disorder triggered by gluten ingestion that damages the small intestine. Although CeD has been associated with a higher risk for cancer, the role of CeD as a risk factor for specific malignancies, such as enteropathy-associated T-cell lymphoma (EATL), remains controversial. Using two-sample Mendelian randomization (2SMR) methods and the summarized results of large genome-wide association studies from public repositories, we addressed the causal relationship between CeD and eight different malignancies. Eleven non-HLA SNPs were selected as instrumental variables (IVs), and causality estimates were obtained using four 2SMR methods: random-effects inverse variance-weighted, weighted median estimation, MR-Egger regression, and MR pleiotropy residual sum and outlier (MR-PRESSO). We identified a significant causal relationship between CeD and mature T/NK cell lymphomas. Under a multivariate Mendelian randomization model, we observed that the causal effect of CeD was not dependent on other known lymphoma risk factors. We found that the most instrumental IV was located in the TAGAP locus, suggesting that aberrant T cell activation might be relevant in the T/NK cell malignization process. Our findings provide new insights into the connection between immune imbalance and the development of severe comorbidities, such as EATL, in patients with CeD.

Immune and spermatogenesis-related loci are involved in the development of extreme patterns of male infertility.

We conducted a genome-wide association study in a large population of infertile men due to unexplained spermatogenic failure (SPGF). More than seven million genetic variants were analysed in 1,274 SPGF cases and 1,951 unaffected controls from two independent European cohorts. Two genomic regions were associated with the most severe histological pattern of SPGF, defined by Sertoli cell-only (SCO) phenotype, namely the MHC class II gene HLA-DRB1 (rs1136759, P = 1.32E-08, OR = 1.80) and an upstream locus of VRK1 (rs115054029, P = 4.24E-08, OR = 3.14), which encodes a protein kinase involved in the regulation of spermatogenesis. The SCO-associated rs1136759 allele (G) determines a serine in the position 13 of the HLA-DRß1 molecule located in the antigen-binding pocket. Overall, our data support the notion of unexplained SPGF as a complex trait influenced by common variation in the genome, with the SCO phenotype likely representing an immune-mediated condition.

The Effect of Body Fat Distribution on Systemic Sclerosis.

Obesity contributes to a chronic proinflammatory state, which is a known risk factor to develop immune-mediated diseases. However, its role in systemic sclerosis (SSc) remains to be elucidated. Therefore, we conducted a two-sample mendelian randomization (2SMR) study to analyze the effect of three body fat distribution parameters in SSc. As instrumental variables, we used the allele effects described for single nucleotide polymorphisms (SNPs) in different genome-wide association studies (GWAS) for SSc, body mass index (BMI), waist-to-hip ratio (WHR) and WHR adjusted for BMI (WHRadjBMI). We performed local (pHESS) and genome-wide (LDSC) genetic correlation analyses between each of the traits and SSc and we applied several Mendelian randomization (MR) methods (i.e., random effects inverse-variance weight, MR-Egger regression, MR pleiotropy residual sum and outlier method and a multivariable model). Our results show no genetic correlation or causal relationship between any of these traits and SSc. Nevertheless, we observed a negative causal association between WHRadjBMI and SSc, which might be due to the effect of gastrointestinal complications suffered by the majority of SSc patients. In conclusion, reverse causality might be an especially difficult confounding factor to define the effect of obesity in the onset of SSc.

Recent advances in elucidating the genetic basis of systemic sclerosis.

PURPOSE OF REVIEW: Systemic sclerosis (SSc) is a complex autoimmune disorder that affects the connective tissue and causes severe vascular damage and fibrosis of the skin and internal organs. There are recent advances in the field that apply novel methods to high throughput genotype information of thousands of patients with SSc and provide promising results towards the use of genomic data to help SSc diagnosis and clinical care. RECENT FINDINGS: This review addresses the development of the first SSc genomic risk score, which can contribute to differentiating SSc patients from healthy controls and other immune-mediated diseases. Moreover, we explore the implementation of data mining strategies on the results of genome-wide association studies to highlight subtype-specific HLA class II associations and a strong association of the HLA class I locus with SSc for the first time. Finally, the combination of genomic data with transcriptomics informed drug repurposing and genetic association studies in well characterized SSc patient cohorts identified markers of severe complications of the disease. SUMMARY: Early diagnosis and clinical management of SSc and SSc-related complications are still challenging for rheumatologists. The development of predictive models and tools using genotype data may help to finally deliver personalized clinical care and treatment for patients with SSc in the near future.

Common genetic variation in KATNAL1 non-coding regions is involved in the susceptibility to severe phenotypes of male infertility.

BACKGROUND: Previous studies in animal models evidenced that genetic mutations of KATNAL1, resulting in dysfunction of its encoded protein, lead to male infertility through disruption of microtubule remodelling and premature germ cell exfoliation. Subsequent studies in humans also suggested a possible role of KATNAL1 single-nucleotide polymorphisms in the development of male infertility as a consequence of severe spermatogenic failure. OBJECTIVES: The main objective of the present study is to evaluate the effect of the common genetic variation of KATNAL1 in a large and phenotypically well-characterised cohort of infertile men because of severe spermatogenic failure. MATERIALS AND METHODS: A total of 715 infertile men because of severe spermatogenic failure, including 210 severe oligospermia and 505 non-obstructive azoospermia patients, as well as 1058 unaffected controls were genotyped for three KATNAL1 single-nucleotide polymorphism taggers (rs2077011, rs7338931 and rs2149971). Case-control association analyses by logistic regression assuming different models and in silico functional characterisation of risk variants were conducted. RESULTS: Genetic associations were observed between the three analysed taggers and different severe spermatogenic failure groups. However, in all cases, the haplotype model (rs2077011*C | rs7338931*T | rs2149971*A) better explained the observed associations than the three risk alleles independently. This haplotype was associated with non-obstructive azoospermia (adjusted p = 4.96E-02, odds ratio = 2.97), Sertoli-cell only syndrome (adjusted p = 2.83E-02, odds ratio = 5.16) and testicular sperm extraction unsuccessful outcomes (adjusted p = 8.99E-04, odds ratio = 6.13). The in silico analyses indicated that the effect on severe spermatogenic failure predisposition could be because of an alteration of the KATNAL1 splicing pattern. CONCLUSIONS: Specific allelic combinations of KATNAL1 genetic polymorphisms may confer a risk of developing severe male infertility phenotypes by favouring the overrepresentation of a short non-functional transcript isoform in the testis.

Common Variation in the PIN1 Locus Increases the Genetic Risk to Suffer from Sertoli Cell-Only Syndrome.

We aimed to analyze the role of the common genetic variants located in the PIN1 locus, a relevant prolyl isomerase required to control the proliferation of spermatogonial stem cells and the integrity of the blood-testis barrier, in the genetic risk of developing male infertility due to a severe spermatogenic failure (SPGF). Genotyping was performed using TaqMan genotyping assays for three PIN1 taggers (rs2287839, rs2233678 and rs62105751). The study cohort included 715 males diagnosed with SPGF and classified as suffering from non-obstructive azoospermia (NOA, n = 505) or severe oligospermia (SO, n = 210), and 1058 controls from the Iberian Peninsula. The allelic frequency differences between cases and controls were analyzed by the means of logistic regression models. A subtype specific genetic association with the subset of NOA patients classified as suffering from the Sertoli cell-only (SCO) syndrome was observed with the minor alleles showing strong risk effects for this subset (ORaddrs2287839 = 1.85 (1.17-2.93), ORaddrs2233678 = 1.62 (1.11-2.36), ORaddrs62105751 = 1.43 (1.06-1.93)). The causal variants were predicted to affect the binding of key transcription factors and to produce an altered PIN1 gene expression and isoform balance. In conclusion, common non-coding single-nucleotide polymorphisms located in PIN1 increase the genetic risk to develop SCO.

Genome-wide compound heterozygote analysis highlights DPY19L2 alleles in a non-consanguineous Spanish family with total globozoospermia.

RESEARCH QUESTION: Would the use of genome-wide genotyping be an advantageous strategy to identify the molecular aetiology of two brothers from a non-consanguineous family, clinically diagnosed with total globozoospermia? DESIGN: Two related Spanish globozoospermic patients were studied. Eight first- and second-degree family members were also included in the study. The clinical procedure included anamnesis, physical examination and semen analyses. Acrosome visualization was performed by fluorescein isothiocyanate-Pisum sativum agglutinin labelling and ultrastructural electron microscope sperm analysis. Sperm DNA fragmentation was determined by TUNEL and SCD. Molecular analysis included: the detection of deletion of the DPY19L2 gene by a BPa (break point "a") gap-polymerase chain reaction, and genotyping by using a high-throughput genome-wide genotyping platform and a genotype imputation strategy. RESULTS: The biological characteristics of the two globozoospermic siblings included round-headed spermatozoa without an acrosome; ultrastructural defects in spermatozoa; increased sperm fragmentation and aneuploidies, inability of spermatozoa to activate oocytes (correctable with artificial activation) and good developmental potential of embryos generated by IVF/intracytoplasmic sperm injection. This genetic study focused on a genome-wide compound heterozygote analysis that identified two deleterious rare coding variants in the DPY19L2 gene [rs771726551 (c.431T>A exon 3) and rs147579680 (c.869G>A exon 8)]. CONCLUSION: A genome-wide compound heterozygote analysis strategy should be considered for molecular screening in globozoospermia and other rare congenital diseases, particularly in cases from non-consanguineous families.

Immune disease risk variants regulate gene expression dynamics during CD4+ T cell activation.

During activation, T cells undergo extensive gene expression changes that shape the properties of cells to exert their effector function. Understanding the regulation of this process could help explain how genetic variants predispose to immune diseases. Here, we mapped genetic effects on gene expression (expression quantitative trait loci (eQTLs)) using single-cell transcriptomics. We profiled 655,349 CD4+ T cells, capturing transcriptional states of unstimulated cells and three time points of cell activation in 119 healthy individuals. This identified 38 cell clusters, including transient clusters that were only present at individual time points of activation. We found 6,407 genes whose expression was correlated with genetic variation, of which 2,265 (35%) were dynamically regulated during activation. Furthermore, 127 genes were regulated by variants associated with immune-mediated diseases, with significant enrichment for dynamic effects. Our results emphasize the importance of studying context-specific gene expression regulation and provide insights into the mechanisms underlying genetic susceptibility to immune-mediated diseases.

Immune disease variants modulate gene expression in regulatory CD4+ T cells.

Identifying cellular functions dysregulated by disease-associated variants could implicate novel pathways for drug targeting or modulation in cell therapies. However, follow-up studies can be challenging if disease-relevant cell types are difficult to sample. Variants associated with immune diseases point toward the role of CD4+ regulatory T cells (Treg cells). We mapped genetic regulation (quantitative trait loci [QTL]) of gene expression and chromatin activity in Treg cells, and we identified 133 colocalizing loci with immune disease variants. Colocalizations of immune disease genome-wide association study (GWAS) variants with expression QTLs (eQTLs) controlling the expression of CD28 and STAT5A, involved in Treg cell activation and interleukin-2 (IL-2) signaling, support the contribution of Treg cells to the pathobiology of immune diseases. Finally, we identified seven known drug targets suitable for drug repurposing and suggested 63 targets with drug tractability evidence among the GWAS signals that colocalized with Treg cell QTLs. Our study is the first in-depth characterization of immune disease variant effects on Treg cell gene expression modulation and dysregulation of Treg cell function.

Contribution of TEX15 genetic variants to the risk of developing severe non-obstructive oligozoospermia.

Background: Severe spermatogenic failure (SPGF) represents one of the most relevant causes of male infertility. This pathological condition can lead to extreme abnormalities in the seminal sperm count, such as severe oligozoospermia (SO) or non-obstructive azoospermia (NOA). Most cases of SPGF have an unknown aetiology, and it is known that this idiopathic form of male infertility represents a complex condition. In this study, we aimed to evaluate whether common genetic variation in TEX15, which encodes a key player in spermatogenesis, is involved in the susceptibility to idiopathic SPGF. Materials and Methods: We designed a genetic association study comprising a total of 727 SPGF cases (including 527 NOA and 200 SO) and 1,058 unaffected men from the Iberian Peninsula. Following a tagging strategy, three tag single-nucleotide polymorphisms (SNPs) of TEX15 (rs1362912, rs323342, and rs323346) were selected for genotyping using TaqMan probes. Case-control association tests were then performed by logistic regression models. In silico analyses were also carried out to shed light into the putative functional implications of the studied variants. Results: A significant increase in TEX15-rs1362912 minor allele frequency (MAF) was observed in the group of SO patients (MAF = 0.0842) compared to either the control cohort (MAF = 0.0468, OR = 1.90, p = 7.47E-03) or the NOA group (MAF = 0.0472, OR = 1.83, p = 1.23E-02). The genotype distribution of the SO population was also different from those of both control (p = 1.14E-02) and NOA groups (p = 4.33-02). The analysis of functional annotations of the human genome suggested that the effect of the SO-associated TEX15 variants is likely exerted by alteration of the binding affinity of crucial transcription factors for spermatogenesis. Conclusion: Our results suggest that common variation in TEX15 is involved in the genetic predisposition to SO, thus supporting the notion of idiopathic SPGF as a complex trait.

Effect and in silico characterization of genetic variants associated with severe spermatogenic disorders in a large Iberian cohort.

BACKGROUND: Severe spermatogenic failure (SpF) represents the most extreme manifestation of male infertility, as it decreases drastically the semen quality leading to either severe oligospermia (SO, <5 million spermatozoa/mL semen) or non-obstructive azoospermia (NOA, complete lack of spermatozoa in the ejaculate without obstructive causes). OBJECTIVES: The main objective of the present study is to analyze in the Iberian population the effect of 6 single-nucleotide polymorphisms (SNPs) previously associated with NOA in Han Chinese through genome-wide association studies (GWAS) and to establish their possible functional relevance in the development of specific SpF patterns. MATERIALS AND METHODS: We genotyped 674 Iberian infertile men (including 480 NOA and 194 SO patients) and 1058 matched unaffected controls for the GWAS-associated variants PRMT6-rs12097821, PEX10-rs2477686, CDC42BPA-rs3000811, IL17A-rs13206743, ABLIM1-rs7099208, and SOX5-rs10842262. Their association with SpF, SO, NOA, and different NOA phenotypes was evaluated by logistic regression models, and their functional relevance was defined by comprehensive interrogation of public resources. RESULTS: ABLIM1-rs7099208 was associated with SpF under both additive (OR = 0.86, p = 0.036) and dominant models (OR = 0.78, p = 0.026). The CDC42BPA-rs3000811 minor allele frequency was significantly increased in the subgroup of NOA patients showing maturation arrest (MA) of germ cells compared to the remaining NOA cases under the recessive model (OR = 4.45, p = 0.044). The PEX10-rs2477686 SNP was associated with a negative testicular sperm extraction (TESE) outcome under the additive model (OR = 1.32, p = 0.034). The analysis of functional annotations suggested that these variants affect the testis-specific expression of nearby genes and that lincRNA may play a role in SpF. CONCLUSIONS: Our data support the association of three previously reported NOA risk variants in Asians (ABLIM1-rs7099208, CDC42BPA-rs3000811, and PEX10-rs2477686) with different manifestations of SpF in Iberians of European descent, likely by influencing gene expression and lincRNA deregulation.

Genomic Risk Score impact on susceptibility to systemic sclerosis.

OBJECTIVES: Genomic Risk Scores (GRS) successfully demonstrated the ability of genetics to identify those individuals at high risk for complex traits including immune-mediated inflammatory diseases (IMIDs). We aimed to test the performance of GRS in the prediction of risk for systemic sclerosis (SSc) for the first time. METHODS: Allelic effects were obtained from the largest SSc Genome-Wide Association Study (GWAS) to date (9 095 SSc and 17 584 healthy controls with European ancestry). The best-fitting GRS was identified under the additive model in an independent cohort that comprised 400 patients with SSc and 571 controls. Additionally, GRS for clinical subtypes (limited cutaneous SSc and diffuse cutaneous SSc) and serological subtypes (anti-topoisomerase positive (ATA+) and anti-centromere positive (ACA+)) were generated. We combined the estimated GRS with demographic and immunological parameters in a multivariate generalised linear model. RESULTS: The best-fitting SSc GRS included 33 single nucleotide polymorphisms (SNPs) and discriminated between patients with SSc and controls (area under the receiver operating characteristic (ROC) curve (AUC)=0.673). Moreover, the GRS differentiated between SSc and other IMIDs, such as rheumatoid arthritis and Sjögren's syndrome. Finally, the combination of GRS with age and immune cell counts significantly increased the performance of the model (AUC=0.787). While the SSc GRS was not able to discriminate between ATA+ and ACA+ patients (AUC<0.5), the serological subtype GRS, which was based on the allelic effects observed for the comparison between ACA+ and ATA+ patients, reached an AUC=0.693. CONCLUSIONS: GRS was successfully implemented in SSc. The model discriminated between patients with SSc and controls or other IMIDs, confirming the potential of GRS to support early and differential diagnosis for SSc.

A distal enhancer at risk locus 11q13.5 promotes suppression of colitis by Treg cells.

Genetic variations underlying susceptibility to complex autoimmune and allergic diseases are concentrated within noncoding regulatory elements termed enhancers1. The functions of a large majority of disease-associated enhancers are unknown, in part owing to their distance from the genes they regulate, a lack of understanding of the cell types in which they operate, and our inability to recapitulate the biology of immune diseases in vitro. Here, using shared synteny to guide loss-of-function analysis of homologues of human enhancers in mice, we show that the prominent autoimmune and allergic disease risk locus at chromosome 11q13.52-7 contains a distal enhancer that is functional in CD4+ regulatory T (Treg) cells and required for Treg-mediated suppression of colitis. The enhancer recruits the transcription factors STAT5 and NF-κB to mediate signal-driven expression of Lrrc32, which encodes the protein glycoprotein A repetitions predominant (GARP). Whereas disruption of the Lrrc32 gene results in early lethality, mice lacking the enhancer are viable but lack GARP expression in Foxp3+ Treg cells, which are unable to control colitis in a cell-transfer model of the disease. In human Treg cells, the enhancer forms conformational interactions with the promoter of LRRC32 and enhancer risk variants are associated with reduced histone acetylation and GARP expression. Finally, functional fine-mapping of 11q13.5 using CRISPR-activation (CRISPRa) identifies a CRISPRa-responsive element in the vicinity of risk variant rs11236797 capable of driving GARP expression. These findings provide a mechanistic basis for association of the 11q13.5 risk locus with immune-mediated diseases and identify GARP as a potential target in their therapy.

Evaluation of Male Fertility-Associated Loci in a European Population of Patients with Severe Spermatogenic Impairment.

Infertility is a growing concern in developed societies. Two extreme phenotypes of male infertility are non-obstructive azoospermia (NOA) and severe oligospermia (SO), which are characterized by severe spermatogenic failure (SpF). We designed a genetic association study comprising 725 Iberian infertile men as a consequence of SpF and 1058 unaffected controls to evaluate whether five single-nucleotide polymorphisms (SNPs), previously associated with reduced fertility in Hutterites, are also involved in the genetic susceptibility to idiopathic SpF and specific clinical entities. A significant difference in the allele frequencies of USP8-rs7174015 was observed under the recessive model between the NOA group and both the control group (p = 0.0226, OR = 1.33) and the SO group (p = 0.0048, OR = 1.78). Other genetic associations for EPSTI1-rs12870438 and PSAT1-rs7867029 with SO and between TUSC1-rs10966811 and testicular sperm extraction (TESE) success in the context of NOA were observed. In silico analysis of functional annotations demonstrated cis-eQTL effects of such SNPs likely due to the modification of binding motif sites for relevant transcription factors of the spermatogenic process. The findings reported here shed light on the molecular mechanisms leading to severe phenotypes of idiopathic male infertility, and may help to better understand the contribution of the common genetic variation to the development of these conditions.

GWAS for systemic sclerosis identifies multiple risk loci and highlights fibrotic and vasculopathy pathways.

Systemic sclerosis (SSc) is an autoimmune disease that shows one of the highest mortality rates among rheumatic diseases. We perform a large genome-wide association study (GWAS), and meta-analysis with previous GWASs, in 26,679 individuals and identify 27 independent genome-wide associated signals, including 13 new risk loci. The novel associations nearly double the number of genome-wide hits reported for SSc thus far. We define 95% credible sets of less than 5 likely causal variants in 12 loci. Additionally, we identify specific SSc subtype-associated signals. Functional analysis of high-priority variants shows the potential function of SSc signals, with the identification of 43 robust target genes through HiChIP. Our results point towards molecular pathways potentially involved in vasculopathy and fibrosis, two main hallmarks in SSc, and highlight the spectrum of critical cell types for the disease. This work supports a better understanding of the genetic basis of SSc and provides directions for future functional experiments.

Chromatin activity at GWAS loci identifies T cell states driving complex immune diseases.

Immune-disease-associated variants are enriched in active chromatin regions of T cells and macrophages. However, whether these variants function in specific cell states is unknown. Here we stimulated T cells and macrophages in the presence of 13 cytokines and profiled active and open chromatin regions. T cell activation induced major chromatin remodeling, while the presence of cytokines fine-tuned the magnitude of changes. We developed a statistical method that accounts for subtle changes in the chromatin landscape to identify SNP enrichment across cell states. Our results point towards the role of immune-disease-associated variants in early rather than late activation of memory CD4+ T cells, with modest differences across cytokines. Furthermore, variants associated with inflammatory bowel disease are enriched in type 1 T helper (TH1) cells, whereas variants associated with Alzheimer's disease are enriched in different macrophage cell states. Our results represent an in-depth analysis of immune-disease-associated variants across a comprehensive panel of activation states of T cells and macrophages.